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Combined analyses of recent cosmological data are showing interesting hints for the presence of 
an extra relativistic component, coined Dark Radiation. Here we perform a new search for Dark 
Radiation, parametrizing it with an effective number of relativistic degrees of freedom parameter, 
Nag. We show that the cosmological data we considered are clearly suggesting the presence for 
an extra relativistic component with N e s = 4.08io'e;g a t 95% c.L. Performing an analysis on Dark 
Radiation sound speed c e a and viscosity c v j s parameters, we found Ce ff = 0.312 ± 0.026 and c? is = 
0. 29tg' 16 a ^ 95% c.L, consistent with the expectations of a relativistic free streaming component 
(Ce ff =cJ is = l/3). Assuming the presence of 3 relativistic neutrinos we constrain the extra relativistic 
component with N„ = I.IOIq ™ and c^ ff = 0.24to'i-j at 95% c.L while c^ is results as unconstrained. 
Assuming a massive neutrino component we obtain further indications for Dark Radiation with 
N? = 1.12i°;^ at 95% c.L . 
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I. INTRODUCTION 

Since almost a decade, observations from Cosmic Mi- 
crowave Background (CMB hereafter) satellite, balloon- 
borne and ground based experiments ([3, 0, 0], 01), 
galaxy rcdshift surveys Q and luminosity distance mea- 
surements, are fully confirming the theoretical predic- 
tions of the standard ACDM cosmological model. This 
not only permits to place stringent constraints on the 
parameters of the model but can be fruitfully used to 
constrain non standard physics at the fundamental level, 
such as classes of elementary particle models predicting 
a different radiation content in the Universe. 

One of the major theoretical predictions of the stan- 
dard scenario is the existence of a relativistic energy com- 
ponent ( see e.g. 0), beside CMB photons, with a cur- 
rent energy density given by : 
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where p 1 is the energy density of the CMB photons 
background at temperature T 7 = 2. 728 A" and N c g is in 
principle a free parameter, defined as the effective num- 
ber of relativistic degrees of freedom. Assuming standard 
electrowcak interactions, three active massless neutrinos 
and including the (small) effect of neutrino flavour oscil- 
lations the expected value is 7V e ff = 3.046 with a devia- 
tion from N e ff = 3 that takes into account effects from 
the non-instantaneous neutrino decoupling from the pri- 
mordial photon-baryon plasma (see e.g. 0). 

In recent years, thanks to the continuous experimental 
advancements, the value of N e g has been increasingly 
constrained from cosmology (fl, @, 0, [H, 0, Ul, 
0, 0, HI, 0, 0, 0), ^ing out N eS = at high 
significance. 

However, especially after the new ACT @ and SPT 
[J] CMB results, the data seem to suggest values higher 
than the "standard" one, with N e g ~ 4 — 5 (see e.g. fl3| |. 
[HL 0> 0: 0) m tension with the expected standard 
value at about two standard deviations. 



The number of relativistic degrees of freedom obviously 
depends on the decoupling process of the neutrino back- 
ground from the primordial plasma. However, a value of 
N c fi = 4 is difficult to explain in the three neutrino frame- 
work since non-standard neutrino decoupling is expected 
to maximally increase this value up to N c g ~ 3.12 (see 
e.g. [13l). A possible explanation could be the existence 
of a fourth (or fifth) sterile neutrino. The hypothesis 
of extra neutrino flavour is interesting since recent re- 
sults from short-baseline neutrino oscillation data from 
LSND [13 and MiniBooNE [H[ experiments are con- 
sistent with a pos sible fourth (or fifth) sterile neutrino 
specie (see [l3l Il4{ and references therein). Moreover, a 
larger value for N c g ~ 4 could arise from a completely 
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different physics, related to axions (see e.g. 
ity waves C I24I0 . decaying particles (see e.g. 
dimensions |26l. [27| and dark energy (see e.g. 
references therein). 

As a matter of fact, any physical mechanism able to 
produce extra "dark" radiation produces the same ef- 
fects on the background expansion of additional neutri- 
nos, yielding a larger value for N D f[ from observations. 

Since there is a large number of models that could en- 
hance TVeff it is clearly important to investigate the pos- 
sible ways to discriminate among them. If Dark Radia- 
tion is made of relativistic particles as sterile neutrinos it 
should behave as neutrinos also from the point of view of 
perturbation theory, i.e. if we consider the set of equa- 
tions that describes perturbations in massless neutrino 
(following the definition presented in [29|): 
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it should have an effective sound speed c ff and a viscos- 
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ity speed c v i s such that c* s = c^ is = 1/3. Free streaming 
of relativistic neutrino will indeed produce anisotropics 
in the neutrino background yielding a value of c^ is = 1/3 
while a smaller value would indicate possible non stan- 
dard interactions (see e.g. Ho[). A value of c V i S different 
from zero, as expected in the standard scenario, has been 
detected in [3l| and confirmed in subsequent papers [32| ■ 
More recently, the analysis of [l!| confirmed the presence 
of anisotropies from current cosmological data but also 
suggested the presence of a lower value for the effective 
sound speed with Cg ff = 1/3 ruled out at more than two 
standard deviations. 

Given the current situation and the experimental hints 
for iV ff ~ 4 is therefore timely to perform a new anal- 
ysis for N e g (and the perturbation parameters c^ s and 
cj is ) with the most recent cosmological data. This is the 
kind of analysis we perform in this paper, organizing our 
work as follows: in Sec. II we describe the data and the 
data analysis method adopted. We present our results 
in the first two subsections of Sec. Ill, depending on 
two adopted different parametrizations for the Dark Ra- 
diation. Moreover a model independent analysis is also 
discussed in the last subsection of Sec. III. Finally we 
conclude in Sec. IV. 



in [||, where the Y p parameter is varied assuming Big 
Bang Nucleosynthesis for each values of N e g and 51b in 
the chain. Since the cosmological epoch and the energy 
scales probed by BBN are dramatically different from the 
ones probed by CMB and large scale structure we pre- 
fer to do not assume standard BBN in our analysis and 
to leave the primordial Helium abundance as fixed to a 
value consistent with current observations. 

We account for foregrounds contributions including 
three extra amplitudes: the SZ amplitude A$z, the am- 
plitude of clustered point sources Ac, and the amplitude 
of Poisson distributed point sources Ap. We marginalize 
the contribution from point sources only for the ACT and 
SPT data, based on the templates provided by [||. We 
quote only one joint amplitude parameter for each com- 
ponent (clustered and Poisson distributed). Instead, the 
SZ amplitude is obtained fitting the WMAP data with 
the WMAP own template, while for SPT and ACT it is 
calculated using the HH SZ template at 148 GHz. A gam, 
this is different from the analysis performed in Q where 
no SZ contribution was considered for the WMAP data. 



III. RESULTS 



II. ANALYSIS METHOD 

We perform a COSMOMC [33| analysis combining the 
following CMB datasets: WMAP 7 Q, ACBAR Q, ACT 
0, and SPT Q, and we analyze datasets using out to 
( max = 3000. We also include information on dark mat- 
ter clustering from the galaxy power spectrum extracted 
from the SDSS-DR7 luminous red galaxy sample Fi- 
nally, we impose a prior on the Hubble parameter based 
on the last Hubble Space Telescope observations [34|. 

The analysis method we adopt is based on the publicly 
available Monte Carlo Markov Chain package cosmomc 
[33j with a convergence diagnostic done through the 
Gelman and Rubin statistic. We sample the following 
six-dimensional standard set of cosmological parameters, 
adopting flat priors on them: the baryon and cold dark 
matter densities 51b and 51 c , the ratio of the sound hori- 
zon to the angular diameter distance at decoupling 8, the 
optical depth to reionization r, the scalar spectral index 
ns, and the overall normalization of the spectrum As at 
k = 0.002 Mpc" 1 . We consider purely adiabatic initial 
conditions and we impose spatial flatness. We vary the 
effective number of relativistic degrees of freedom N e g, 
the effective sound speed c^ ff , and the viscosity param- 
eter c^ is . In some cases, we consider only variations in 
the extra dark radiation component N„ = N e g — 3.046, 
varying the perturbation parameters c V i S and c e ff only for 
this extra component and assuming c^ s = c^ is = 1/3 for 
the standard 3 neutrino component. 

In our analysis we always fix the primordial Helium 
abundance to the observed value Y p = 0.24. This pro- 
cedure is different from the one adopted, for example, 



As stated in the previous section, we perform two dif- 
ferent analyses. In the first analysis we vary the ampli- 
tude of the whole relativistic contribution changing iV e ff 
and the corresponding perturbation parameters c^ is and 
c^ ff . In the second analysis we assume the existence of 
a standard neutrino background and vary only the extra 
component = N e g — 3.046 considering only in this 
extra component the variations in c^ is and c^ ff . 

A. Varying the number of relativistic degrees of 
freedom N e g . 

In TableUwe report the constraints on the cosmological 
parameters varying N e g with and without variations in 
perturbation theory. We consider two cases: first we run 
our analysis fixing the perturbation parameters to the 
standard values, i.e. c^ s = cj is = 1/3, then we let those 
parameters to vary freely. 

As we can see from the results in the left column of 
Table H the WMAP7+ACT+SPT+DR7+H0 analysis is 
clearly suggesting the presence for Dark Radiation with 
N e g = 4.08lo'gg & t c.L. When considering variations 
in the perturbation parameters (right column) the con- 
straint is somewhat shifted towards smaller values with 
A c ff = 3.89^q'7q. The constraint on the sound speed, 
c^ ff = 0.312 ± 0.026 is fully consistent with the expec- 
tations of a free streaming component. Anisotropies in 
the neutrino background are detected at high statisti- 
cal significance with c^ is = 0.29^°'^ improving previous 
constraints presented in [3l| . 

It is interesting to consider the possible degeneracies 
between N c g and other "indirect" (i.e. not considered 
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FIG. 1. 68% and 95% c.l. constraints for the degeneracy between N e g and the Hubble constant Ho, the age of the universe to, 
and the amplitude of mass fluctuations erg. 
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0.02229 ± 0.00038 
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TABLE I. MCMC estimation of the cosmological parameters 
assuming JV c ff relativistic neutrinos. Upper bounds at 95% c.l. 
are reported for foregrounds parameters. We quote the one- 
dimensional marginalized 68% and 95% c.l. for the neutrino 
parameters. 



as primary parameters in MCMC runs) model param- 
eters. In Figure [T] we therefore plot the 2D likelihood 
constraints on N e g versus the Hubble constant Ho, the 
age of the universe to and the amplitude of r.m.s. mass 
fluctuations on spheres of 8Mpc/i _1 , eg. 

As we can see from the three panels in the figure, there 
is a clear degeneracy between N e g and those three param- 
eters. Namely, an extra radiation component will bring 
the cosmological constraints (respect to the standard 3 
neutrino case) to higher values of the Hubble constant 
and of a$ and to lower values of the age of the universe 
t . These degeneracies have been already discussed in 
the literature (see e.g. (36[) and could be useful to esti- 
mate the effect of additional datasets on our result. The 
3% determination of the Hubble constant from the anal- 
ysis of p~7j | plays a key role in our analysis in shifting 
the constraints towards larger values of iV c ff. If future 



analyses will point towards lower values of the Hubble 
constant, this will make the standard 3 neutrino case 
more consistent with observations. If future observations 
will point towards values of the age of the universe signif- 
icantly larger than 13 Gyrs, this will be against an extra 
dark radiation component, since it prefers to ~ 12.5Gyrs. 
Clearly, adding cluster mass function data as presented 
in (37$ and that points towards lower values of as ren- 
ders the standard N e s = 3.046 case more consistent with 
observations. A future and precise determination of erg 
from clusters or Lyman-a surveys could be crucial in rul- 
ing out dark radiation. 



B. Varying only the excess in the relativistic 
component Nj, and assuming 3 standard neutrinos. 

In Table ITTI we report the constraints considering only 
an excess N„ in the number of relativistic degrees of 
freedom over a standard 3 neutrinos background. 

As we can see for the results in the table, the evidence 
for an extra background is solid with = 1.46^0 74 
at 95% c.l. when only variations in the c^ is component 
are considered, while the constraint is N„ = l.lO+o 72 
when also variations in Cg ff are considered. Again, the 
data provide a good determination for c^ ff with Cg ff = 
0.24^q'53 at 95% c.l., in marginal agreement at about 2a 
with the standard c^ ff = 1/3 value. This lower value of 
c^ ff , also found in [19J, could hint for a dark radiation 
component with a varying equation of state, ruling out a 
a masslcss sterile neutrino. It will be certainly interesting 
to investigate if this signal remains in future analyses. No 
significant constraint is obtained on c^ is . 

In Figure [2] we show the degeneracy between the pa- 
rameters iVjf , Cg ff , and c^ is by plotting the 2D likelihood 
contours between them. As we can see a degeneracy is 
present between c^ ff and Njjj: models with lower values 
of N„ are more compatible with c^ ff = since the effect 
of Cg ff on the CMB spectrum is smaller. No apparent 
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FIG. 2. 68% and 95% c.l. constraints for the degeneracy between neutrinos parameters. Red contours refer to model (A) in 
Table [TTJ while blue contours show model (B). 
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TABLE II. MCMC estimation of the cosmological parameters 
considering an extra component N„ and assuming a standard 
background of 3 relativistic neutrinos. The perturbation pa- 
rameters refer to the extra component. Both 68% and 95% 
confidence levels for the neutrino parameters are reported. 
Upper bounds are at 95% c.l. . 



degeneracy is present between c^ is and the remaining pa- 
rameters since c^ is is weakly constrained by current data. 

Since oscillation experiments have clearly established 
that neutrino are massive, it is interesting to perform a 
similar analysis but letting the 3 neutrino standard back- 
ground with c^ ff = c^ is = 1/3 to be massive, and varying 
the parameter Em^ that consider the sum of masses of 
the 3 active neutrinos. The extra dark radiation compo- 
nent is assumed massless and we treat the perturbations 
in it as in the previous sections. In Table Hill we report 
the results of this analysis. 

As we can see, when masses in the active neutrinos 
are considered, there is a slightly stronger evidence for 
the extra background with N„ = 1.12* q'^®. This is can 
be explained by the degeneracy present between ^2 Tri v 
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TABLE III. MCMC estimation of the cosmological param- 
eters considering N v — 3.04 massive neutrinos. Values and 
68% - 95% errors for the neutrino parameters are reported. 
Upper bounds are at 95% c.l. . 



and N„ , well known in the literature (see e.g. [13j) 
and clearly shown in Figure where we report the 2D 
marginalized contours in the plane ^2 m f — ■ 



C. Profile likelihood analysis 

Recently, in [38j , a model- independent analysis for the 
extra relativistic degrees of freedom in cosmological data 
has been performed claiming no statistically significant 
evidence for it. This simple analysis consists in extracting 
the maximum likelihood value L as a function of N e s over 
the parameter space sampled in the chains, with a bin 
width of 0.5 and constructing a profile likelihood ratio by 
considering ln(LN oll / L max ) as a function of N s; where 
L m ax is the maximum likelihood in the entire chains. 

Here we perform a similar analysis, using however a 
smaller bin width of 0.05 and considering the case where 
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0.5 1.0 1.5 2.0 2.5 3.0 

N S 
v 

FIG. 3. Degeneracy in the plane ^ m v — N„ at 68% and 95% 
c.l. . 




FIG. 4. Maximum Likelihood ratio I/jv off / ' L max for N e g. The 
dashed lines represent the 68% and 95% c.l. for a Gaussian 
likelihood (L Nott /L max = 0.6065 and L No{{ /L max = 0.135) 
respectively. 



the whole number of relativistic degrees of freedom N e g 
is varied while c 2 is = c 2 ff = 1/3. The resulting likelihood 
ratio Ljv eff / L max , plotted in Figure |4j clearly indicates 
a preference for a dark radiation component finding that 
the best fit model has N eS = 3.88 with a Ax 2 = 14.56 



respect to the best fit model with N c g = 3.046. 

We should however point out that the ratio Ln b{{ / L max 
presented in Figure [4] is rather noisy. Bayesian methods 
such as MCMC are indeed known to be inaccurate for this 
purpose (see for example the discussion in [4l[). Other 
methods more appropriate for a frequentist analysis have 
been presented, for example, in [42^ . 

IV. CONCLUSIONS 

In this paper we performed a new search for Dark Ra- 
diation, parametrizing it with an effective number of rela- 
tivistic degrees of freedom N c g. We have shown that the 
cosmological data we considered are clearly suggesting 
the presence for an extra dark radiation component with 
N fi — 4.08^ 68 & t 95% c.l. . Performing an analysis on 
its effective sound speed c e g and viscosity c v - ls parame- 
ters, we found c 2 cff = 0.312 ± 0.026 and c 2 is = 0.29±g;^ 
at 95% c.l., consistent with the expectations of a rela- 
tivistic free streaming component (c 2 ff =c 2 is =l/3). As- 
suming the presence of 3 standard relativistic neutrinos 
we constrain the extra dark radiation component with 

N$ = 1.101:8;™ and c off = °- 24 ±ai3 at 95% C - L whilc 
c 2 is is practically unconstrained. Assuming a mass in the 
3 neutrino component we obtain further indications for 
the dark radiation component with = 1.12lto74 a t 
95% c.l. . From these results we conclude that Dark 
Radiation currently represents one of the most relevant 
anomaly for the A-CDM scenario. 

When comparison is possible, our results are in good 
agreement with the most recent analysis presented in [l9[ 
that uses a different choice of datasets (for example, we 
don't consider matter fluctuations data from Lyman-a as 
in [13]) and an independent analysis method. 

Dark Radiation will be severely constrained in the 
very near future by the Planck satellite data, where a 
precision on N e g of about A7V ff ~ 0.2 is expected (see 
e.g. [Hi and 0) only from CMB data. 
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